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Abstract 

This paper is a revised version of our recent publication Faber et al, Phys. Rev. 
D62 (2000) 025019, hep-th /990704£ . The main revision concerns the expansion 
into group characters that we have used for the evaluation of path integrals over 
gauge degrees of freedom. In the present paper we apply the expansion recom- 

OO ' 
O 



mended by Diakonov and Petrov in hep-lat / 0008004 . Our former expansion was 
£NJ ■ approximate and in the region of the particular values of parameters violated the 

completeness condition by 1.4%. We show that by using the expansion into charac- 



ters recommended by Diakonov and Petrov in hep-fat /0008004 our previous results 
are retained and the path integral over gauge degrees of freedom for Wilson loops 
derived by Diakonov and Petrov (Phys. Lett. B224 (1989) 131 and, correspond- 



ed ,| ingly, |hep-lat / 0008004 ) by using a special regularization is erroneous and predicts 



zero value for the Wilson loop. This property is obtained by direct evaluation of 
path integrals for Wilson loops defined for pure SU(2) gauge fields and Z(2) center 
vortices with spatial azimuthal symmetry. Further we show that both derivations 
^ ■ given by Diakonov and Petrov for their regularized path integral, if done correctly, 

predict also zero value for Wilson loops. Therefore, the application of their path in- 
tegral representation of Wilson loops cannot give a new way to check confinement in 
lattice as has been declared by Diakonov and Petrov (Phys. Lett. B242 (1990) 425 
and hep-lat / 0008004| ) . Our statement pointing out that none non-Abelian Stokes 



theorem can exist for Wilson loops except the old-fashioned one derived by means 
of the path-ordering procedure is retained. It is based on well-defined properties 
of group characters and is not related to whatever explicit method of evaluation of 
path integrals is applied. Comments on the paper |hep-lat/0008004| by Diakonov 
and Petrov are given. Some misprints in our paper Phys. Rev. D62 (2000) 025019 
are corrected. 
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1 Introduction 



The hypothesis of quark confinement, bridging the hypothesis of the existence of 
quarks and the failure of the detection of quarks as isolated objects, is a challenge for 
QCD. As a criterion of colour confinement in QCD, Wilson [1] suggested to consider the 
average value of an operator 

W{C) = ^trVce'9 § C dx,A,{x) = _L t rtf(C ra ), (1.1) 

defined on an closed loop C, where A^(x) = t a A®(x) is a gauge field, t a (a = 1, . . . , N 2 — 1) 
are the generators of the SU(N) gauge group in fundamental representation normalized 
by the condition tr (t a t b ) = S ab /2, g is the gauge coupling constant and Vc is the operator 
ordering colour matrices along the path C. The trace in Eq.( |l.l|) is computed over colour 
indices. The operator 

U{C yx ) = VcJ 9 /c - dz » A ^ z \ (1.2) 

makes a parallel transport along the path C yx from x to y. For Wilson loops the contour C 
defines a closed path C xx . For determinations of the parallel transport operator U(C yx ) the 
action of the path-ordering operator Vc is defined by the following limiting procedure 
[2] 

U{C yx ) = V C J 9 $c «* dZ " Mz) = hm ft tf(CWx) = 

= )imU{C yXn _ x ) . . . U(C X2X1 ) U(C X1X ) = Jim f[ c ^(** " ' A ( x ^) , (1.3) 

k=l 

where C XkXk _ 1 is an infinitesimal segment of the path C yx with xq = x and x n = y. The 
parallel transport operator U{C XkXk _^) for an infinitesimal segment C XkXk _ x is defined by 

[2]: 

U{C XkXk _ x ) = e 9 ^ c vk-i dz ^ A ^ = e i 9 fa - Xk-t) • A(x k -.!) ^ ( L4 ) 

In accordance with the definition of the path-ordering procedure ( |1.3| ) the parallel trans- 
port operator U{C yx ) has the property 

U(C yx ) = U(C yXl )U(C XlX ) } (1.5) 

where x\ belongs to the path C yx . Under gauge transformations with a gauge function 

A^z) - A n (z) = Q(z)A ft (z)^{z) + - d„n(z)tf(z), (1.6) 

W 

the operator U(C yx ) has a very simple transformation law 

U(C yx ) - U n (C yx ) = Q(y) U(C yx ) Q\x). (1.7) 
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We would like to stress that this equation is valid even if the gauge functions Q(x) and 
fl(y) differ significantly for adjacent points x and y. 

As has been postulated by Wilson [1] the average value of the Wilson loop < W(C) > 
in the confinement regime should show area-law falloff [1] 

< W{C) >^e~ aA , (1.8) 

where a and A are the string tension and the minimal area of the loop, respectively. As 
usually the minimal area is a rectangle of size L x T. In this case the exponent a A can 
be represented in the equivalent form o A = V(L) T, where V(L) = oL is the interquark 
potential and L is the relative distance between quark and anti-quark. 

The paper is organized as follows. In Sect. 2 we discuss the path integral representation 
for Wilson loops by using well-known properties of group characters. The discretized form 
of this path integral is naturally provided by properties of group characters and does 
not need any artificial regularization. We derive a closed expression for Wilson loops in 
irreducible representation j of SU (2). In Sect. 3 we extend the path integral representation 
to the gauge group SU(N). As an example, we give an explicit representation for Wilson 
loops in the fundamental representation of ££7(3). In Sects. 4 and 5 we evaluate the path 
integral for Wilson loops, suggested in Ref.[3], for two specific gauge field configurations 
(i) a pure gauge field in the fundamental representation of SU(2) and (ii) Z(2) center 
vortices with spatial azimuthal symmetry, respectively. We show that this path integral 
representation fails to describe the original Wilson loop for both cases. In Sect. 6 we 
show that the regularized evolution operator in Ref. [3] representing Wilson loops in the 
form of the path integral over gauge degrees of freedom has been evaluated incorrectly 
by Diakonov and Petrov. The correct value for the evolution operator is zero. This 
result agrees with those obtained in Sects. 4 and 5. In Sect. 7 we criticize the removal 
of the oscillating factor from the evolution operator suggested in Ref. [3] via a shift of 
energy levels of the axial-symmetric top. We show that such a removal is prohibited. It 
leads to a change of symmetry of the starting system from £77(2) to £7(2). Keeping the 
oscillating factor one gets a vanishing value of Wilson loops in agreement with our results 
in Sects. 4, 5 and 6. In the Conclusion we discuss the obtained results. In Appendix we 
give comments on the paper hep-lat / 0008004 by Diakonov and Petrov. 



2 Path integral representation for Wilson loops 

Attempts to derive a path integral representation for Wilson loops (|1 . 1|) , where the 
path ordering operator is replaced by a path integral, have been undertaken in Refs.[3-5]. 
The path integral representations have been derived for Wilson loops in terms of gauge 
degrees of freedom (bosonic variables) [3,4] and fermionic degrees of freedom (Grassmann 
variables) [5]. For the derivation of the quoted path integral representations for Wilson 
loop different mathematical machineries have been used. Below we discuss the derivation 
of the path integral representation for Wilson loops in terms of gauge degrees of freedom 
by using well-known properties of group characters. In this case a discretized form of path 
integrals is naturally provided by the properties of group characters and the completeness 
condition of gauge functions. It coincides with the standard discretization of Feynman 
path integrals [6] and does not need any artificial regularization. 
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We argue that the path integral representation for Wilson loops suggested by Diakonov 
and Petrov in Ref. [3] is erroneous. For the derivation of this path integral representation 
Diakonov and Petrov have used a special regularization drawing an analogy with an axial- 
symmetric top. The moments of inertia of this top are taken finally to zero. As we show 
below this path integral amounts to zero for Wilson loops defined for 577(2). Therefore, 
it is not a surprise that the application of this erroneous path integral representation 
to the evaluation of the average value of Wilson loops has led to the conclusion that 
for large loops the area-law falloff is present for colour charges taken in any irreducible 
representation r of SU(N) [7]. This statement has not been supported by numerical 
simulations within lattice QCD [8]. As has been verified, e.g. in Ref. [8] for SU(3), in the 
confined phase and at large distances, colour charges with non-zero iV-ality have string 
tensions of the corresponding fundamental representation, whereas colour charges with 
zero iV-ality are screened by gluons and cannot form a string at large distances. Hence, 
the results obtained in Ref. [7] cannot give a new way to check confinement in lattice as 
has been declared by Diakonov and Petrov. 

For the derivation of Wilson loops in the form of a path integral over gauge degrees of 
freedom by using well-known properties of group characters it is convenient to represent 
W(C) in terms of characters of irreducible representations of SU(N) [9-11] 

W r (C) = j X [Ur(C xx )}, (2.1) 

where the matrix U r (C xx ) realizes an irreducible and e^-dimensional matrix representation 
r of the group SU(N) with the character x[U r (C xx )} = tr[U r (C xx )]. 

In order to introduce the path integral over gauge degrees of freedom we suggest to 

use 

J DQ rX [U r Ql] xKVr] = j xWrVr], (2.2) 

where the matrices U r and V r belong to the irreducible representation r, and DQ r is 
the Haar measure normalized to unity / DQ r = 1. The orthogonality relation for gauge 
functions Q r reads 

/ ™ r (nt) aibi (n r ) a2b2 = i s aih 5 bia2 . (2.3) 

By using the orthogonality relation it is convenient to represent the Wilson loop in the 
form of the integral 

W r (C) = j J DQ r (x) X [n r (x)U r (C xx )Ql(x)]. (2.4) 

According to Eq. Ql.3D and Eq. Ql.5D the matrix U r {C xx ) can be decomposed in 

U T (C XX ) = lim U r {C XXn _ 1 )U r {C Xn _ 1Xn _ 2 ) . . . U r (C X2Xl )U r (C XlX ). (2-5) 
Substituting Eq.(|2.5p in Eq.( |2.4[ ) and applying {n — l)-times Eq.(|2.2j) we end up with 

W r (C) = — lirn J ... J 7JO r (xi) . . . fi r (x n ) ^^[^(^^(C^-J^^-i)] 

. . . d rX [n r (xi)u r (c xlXn )ni(x n )}. (2.6) 
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Using relations VL r (x k )U r (C XkXk _ 1 )VLl(x k ^i) = ^(C^^J we get 

W r (C) = ± Jim J . . . J DQ r ( Xl ) . . . DQ r (x n ) d rX [U?(C XnXn ^)} . . . d rX [U?(C xlXn )}.(2.7) 

The integrations over Q r (x k ) (k — 1, . . . , n) are well defined. These are standard integra- 
tions on the compact Lie group SU(N). 

We should emphasize that the integrations over Q r (x k ) (k — 1, . . . , n) are not corre- 
lated and should be carried out independently. 

Due to Eq. ( |2.2|) the discretization of Wilson loops given by Eqs. ( |2.6|) and ( |2.7|) repro- 
duces the standard discretization of Feynman path integrals [6] where infinitesimal time 
steps can be described by a classical motion. Therefore, the discretized expression ( [2.7|) 
can be represented formally by 

Wr(C) = 1 / 1] [d r Dn r (x)\ X [U?(C XX )\. (2.8) 

a r J x eC 

Conversely the evaluation of this path integral corresponds to the discretization given 
by Eqs.( |2.6|) and ( |2.7|) . The measure of the integration over Q r (x) is well defined and 
normalized to unity 

J J] DQ r (x) = Jim J DQ r (x n ) J DVL r (x n ^) ...J Dtt r { Xl ) = 1. (2.9) 

Thus, for the determination of the path integral over gauge degrees of freedom ( |2.8| ) we 
do not need to use any regularization, since the discretization given by Eqs. (|2.6| ) and ( |2.7|) 
are well defined. 

We would like to emphasize that Eq. (|2.8|) is a continuum analogy of the lattice version 
of the path integral over gauge degrees of freedom for Wilson loops used in Eq.(2.13) of 
Ref. [11] for the evaluation of the average value of Wilson loops in connection with Z(2) 
center vortices. 

Now let us to proceed to the evaluation of the characters x[U^(C XkXk _ 1 )]. Due to the 
infinitesimality of the segments C XkXk _ 1 we can omit the path ordering operator in the 
definition of U^{C XkXk _^) [2]. This allows us to evaluate the character x[^?(Ca: fc a; fe _i)] with 
[/^(C^^J taken in the form [2] 

t^(CWJ = expig f dXpAfa). (2.10) 



Of course, the relation given by Eq.( 2.10 ) is only defined in the sense of a meanvalue 
over an infinitesimal segment C XkXk _ 1 . Therefore, it can be regarded to some extent as 
a smoothness condition. Unlike the smoothness condition used by Diakonov and Petrov 
[3] Eq. (|2.10|) does not corrupt the Wilson loop represented by the path integral over the 
gauge degrees of freedom. 

The evaluation of the characters of C^(C f a . ka . fc _ 1 ) given by Eq. (|2.10|) runs as follows. 
First let us consider the simplest case, the 577(2) gauge group, where we have r — j — 
0, 1/2, 1, . . . and dj = 2j + 1. The character XpfiC^^)] is equal to [9,10,12] 

x[U?{C XkXk _ x )] = E <m J \Uf(C XkXk _ 1 )\m J >= 
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^2 ( < >„ , <\> ([,,,,, : :,r* _ (211) 



m j=-3 



where rrij is the magnetic colour quantum number, \rrij > and rrij $[C XfeXfe _ 1 ; A ] are the 
eigenstates and eigenvalues of the operator 

*[C XkXk _ l] A n ] = g [ dx.A^x), (2.12) 

i.e. &[C XkXk _ 1 ; A n ] \rrij >= rrij $[C XfcXfc _ 1 ; A n ] \rrij(xk-i) >• The standard procedure for 
the evaluation of the eigenvalues gives Q[C XkXk _ 1 ; A n ] in the form 

§[C XkXk _ 1 -A n ] = g f y/ gia ,[A a ]{x)dx,>dx v , (2.13) 

where the metric tensor can be given formally by the expression 

9llu [A n ](x) = 2tr[A^](x). (2.14) 

In order to find an explicit expression for the metric tensor we should fix a gauge. As an 
example let us take the Fock-Schwinger gauge 

x^(x) = 0. (2.15) 

In this case the gauge field A^x) can be expressed in terms of the field strength tensor 
G^ u (x) as follows 

i 

Afj,(x) = J ds s x a G a ^(xs). (2.16) 
o 

This can be proven by using the obvious relation 

%aGaii{p£) •^a'^a A^(x) x a 0^A a (x) ig\x a A a (x*), yl^(x)] 

d 

= A^(x) + x a — — An(x), (2.17) 
ox a 

valid for the Fock-Schwinger gauge x a A a (x) = 0. Replacing x — > xs we can represent the 
r.h.s. of Eq.flUTD total derivative with respect to s 

d d 

sx a G afl (xs) = A^(xs) + Xa-^-A^xs) = —[sA^xs)]. (2.18) 

Integrating out s G [0, 1] we arrive at Eq. fl2.16Q . 

Using Eq.( |2.16| ) we obtain the metric tensor g^[A n ](x) in the form 

l l 

g^[A n }{x) = 2x aXf3 f f dsds'ss'tT[G^Jxs)G%(xs')} = 





1 1 



2x a x l3 j J dsds ss'ti[VL(xs)G a ^{xs)Vt\xs)Vt(xs' )G (3v {xs' )D)(xs' )}. (2.19) 





For the derivation of Eq.( 2.19f) we define the operator Q[C XkXk _ 1 ; A n ] of Eq.( 2.12j ) following 



the definition of the phase of the parallel transport operator U(C XkXk _ 1 ) given by Eq.(|T 
[2] 

®[C XkXk _ r ] A n ] = g[ dx ll A^{x) = {x k -x k _ 1 ) ll A^{x k _ 1 ) = 

1 

= (x k -x k -i)^J d88a^_ 1 G^ ll (x k - 1 s). (2.20) 



The parameter s is to some extent an order parameter distinguishing the gauge functions 
Q(x k ) and e nteri ng the relation £l(x k )U(C XkXk _ 1 )£tf(x k -i) = U n (C XkXk _A. 

Substituting Eq.( p.ll ) in Eq. (|2.7p we arrive at the expression for Wilson loops defined 
for SU(2) 

WAG) = t lim 

n ' (2j + l) 2 ™^~ 



J DQj (x n ) (2 j + 1) e ig m ^ Ic *i*n V9^W) dx » dx » 

(n) 

3 J 

fDClj^pj + l) ± e i9m ^ X)$G ^ ^9^[A n ]{x)dx^dx u 



(n-l) 



J DQjixj) (2j + 1) E e 97 "^ ^A An K x ) dx^dx^ ^ ^ 



(!) 
3 J 



The magnetic quantum number mf^ (k = 1, . . . ,n) belongs to the infinitesimal segment 
C Xk+lXk , where C Xn+1Xn (",,,,,.. 

In compact form Eq. fl2.2ip can be written as a path integral over gauge functions 

Wj(C) = 1 J J] DSlfc) J2 (2j + 1) e 9§cm ^ ^d^M^dx^dx^ ^ ^ 

' 3 ' x£C {m,j(x)} 

The integrals along the infinitesimal segments C XkXk _ 1 we determine as [2] 



m.j(x) \jg^{A n ](x) dx^dx u = m^x^x) yj ' g^[A n ](x k -i) Ax^ Ax u 



m 



(fe-i) 



g^[A^]{x k ^)Ax^Ax v . (2.23) 



where Ax = x k — x k ~i- 

Comparing the path integral ( |2.22| ) with that suggested in Eq. (23) of Ref . [3] one finds 
rather strong disagreement. First, this concerns the contribution of different states rrij of 
the representation j. In the case of the path integral ( |2.22| ) there is a summation over all 
values of the magnetic colour quantum number rrij, whereas the representation of Ref. [3] 
contains only one term with rrij = j. Second, Ref. [3] claims that in the integrand of 
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their path integral the exponent should depend only on the gauge field projected onto 
the third axis in colour space. However, this is only possible if the gauge functions are 
slowly varying with x, i.e. fl(xk)^ (xk-i) — 1. In this case the parallel transport operator 
U Q {C XkXk _ x ) would read [13] 

U n (C XkXk _ 1 )=expig f dx ll A%(x) = l + tg{x l -x^ 1 )-A n (x i - 1 ) } (2.24) 

Cx k x k-i 

and the evaluation of the character x\Uf{C XkXk _^\ would run as follows 

< mj\[Uf{C x ^ k _A\\mj >= 1 + (/;;.)„,,„,///(./'/,• - x fc _i) ■ [A n (a*_i)]W = 

= l + m^^(x fc -x fe _ 1 )-[^(x fc _ 1 )]( 3 ) = e yJc Wi " jl n ^ n , (2.25) 

where we have used the matrix elements of the generators of SU (2), i.e. (t°j) mjm . = rrij 5 a3 . 
More generally the exponent on the r.h.s. of Eq.( [2.25|) can be written as 



f dx^mjix) [A^(x)] {3) = 2 f dx^m^x) tr[t|A^(x)]. (2.26) 

J Cx k X k _-y JGx*~<-' 



k X k-l 



This gives the path integral representation for Wilson loops defined for SU(2) in the 
following form 



Wj{C) = 1 / UDnfc) J2(2j + l)e 2i ^cd^m j (x)tx[t 3 j A^(x)] 

(2j + 1) J x( - c {rrijix)} 

(2.27) 



The exponent contains the gauge field projected onto the third axis in colour space 



tr[t^A^(x)]. Nevertheless, Eq.( 2.27 ) differs form Eq. (23) of Ref.[3] by a summation over all 



values of the colour magnetic quantum number rrij of the given irreducible representation 
J- 

The repeated application of Eq. (|2.2j ) induces that the integrations over the gauge func- 
tion at Xk are completely independent of the integrations at Xk±\- There is no mechanism 
which leads to gauge functions smoothly varying with Xk {k = 1, . . . n). In this sense the 
situation is opposite to the quantum mechanical path integral. In Quantum Mechanics 
the integration over all paths is restricted by the kinetic term of the Lagrange function. 
In the semiclassical limit h — > due to the kinetic term the fluctuations of all trajectories 
are shrunk to zero around a classical trajectory. However, in the case of the integration 
over gauge functions for the path integral representation of the Wilson loop, there is 
neither a suppression factor nor a semiclassical limit like h — > 0. The key point of the 
application of Eq. (|2.2j ) and, therefore, the path integral representation for Wilson loops is 
that all integrations over Q(xk) (k — 1, . . . ,n) are completely independent and can differ 
substantially even if the points, where the gauge functions fl(xk) and fi(x&_i) are defined, 
are infinitesimally close to each other. 

For the derivation of Eq.(23) of Ref.[3] Diakonov and Petrov have used at an in- 
termediate step a regularization drawing an analogy with an axial-symmetric top with 
moments of inertia I± and Within this regularization the evolution operator repre- 
senting Wilson loops has been replaced by a path integral over dynamical variables of this 
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axial-symmetric top which correspond to gauge degrees of freedom of the non-Abelian 
gauge field. The regularized expression of the evolution operator has been obtained in 
the limit — * 0. The moments of inertia have been used as parameters like h — > 0. 
Unfortunately, as we show in Sect. 6 the limit I±, Jii — > has been evaluated incorrectly. 



3 The SU(N) extension 



The extension of the path integral representation given in Eq. (|2.24|) to SU(N) is rather 



straightforward and reduces to the evaluation of the character of the matrix U^ l (C XkXk _ 1 ) 
in the irreducible representation r of SU(N). The character can be given by [12] 

= Et^^-^m;^, (3.i) 

m r 

where Hg {i — 1, . . . , N — 1) are diagonal d r x d r traceless matrices realizing the representa- 
tion of the Cartan subalgebra, i.e. [Hi, Hgi] = 0, of the generators of the SU (N) [12]. The 
sum runs over all the weights m T = (m rl , . . . ,m r ff-i) of the irreducible representation r 
and 7m r is the multiplicity of the weight m r and J2m r lm r = d r . The components of the 
vector $[C a ; fcXfc _ 1 ; A n ] are defined by 

*i[C* kn - 1 ;A a \=g f wHx)], (3.2) 

where we have introduced the notation u(x) = t a uo a {x) = dz ■ A n (x). The functions 
ipi [co(x)] are proportional to the roots of the equation det uj{x) — A =0. 

The path integral representation of Wilson loops defined for the irreducible represen- 
tation r of SU (N) reads 

Wr(C) = i / II DQ r (x) £ d rlrhr(x) e^ § c ™r{x)-(p[u{x)]_ (3 _ 3) 

a r J xeC {m r (x)} 

Let us consider in more details the path integral representation of Wilson loops defined 
for the fundamental representation 3 of SU(3). The character X3W3 (C)] is defined as 

X3 [U°(C)] = tv(e iH ^ C ' A ^ +iH2MC;A a \ ) = -i^[C;A n ]/3 

+e i$i[C; A n ]/2V3 e m[C; A n ]/6 + e -i$ a [C; A^v^ e z$ 2 [C; A n ]/6 ) (3 4) 

where ifi = t 3 /\/3 and i? 2 = t 8 /y/3 [12]. For the representation 3 of SU(3) the equation 
det [uj — A] = takes the form 

A 3 - A^W(x) -detcu(x) = 0. (3.5) 
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The roots of Eq.( |3.5|) read 



A (1) 



-— j= \Jtiuo 2 {x) cos y- arccos^ 



2det 



1 + 12 



t a ti{t a uj 2 {x)) 
tr u 2 (x) 



- Jtr(jj 2 (x) sin -arccos,. 



2det 



1 + 12 



t a ti{t a u 2 (x)) 
tr uj 2 (x) 



A (2) 



_ ,/trci; 2 (x) cos f -arccos,, 

s/E v w V 3 A 



2det 



1 + 12 



t a tr(t a uj 2 (x)) 
tr u; 2 (a;) 



— = xlinJHx) sin f - arccos 
y/2 V K J V 3 



\ 



2det 



1 + 12 



t a ti{t a uj 2 (x)) 
trui 2 (x) 



A (3) 



v 'tru; 2 (x) cos I - arccos 



\ 



2det 



1 + 12 



t a tr(t a cj 2 (x)) 
tr u 2 (x) 



In terms of the roots (i — 1, 2, 3) the phases $1,2 [C; v4 Q ] are defined as 



$i[C7;A n ] = -(/V6/ /tr" 



W x sin - arccos 
c • v ; I 3 



$ 2 [C; A n ] = -gVft j> ^ti uo 2 (x) cos Qarccos^ 2det 



1 + 12 



(3.6) 




where trcj 2 (x) = | dx^dx u . Thus, in the fundamental representation 

path integral representation for Wilson loops reads 

W = -/nKW x 3] ( e »i[C;An]/2^ e i$ 2 [C;A°]/6 
9 •/ ,- v . ^ 



+e" 



i^[C; A n ]/2V3 p i<$> 2 [C; A n ]/6 + e -z$ 2 [C; A 



(3i 



where the phases $1,2 [C; A n ] are given by Eq. (|3~7|) . 



4 Wilson loop for pure gauge field 

As has been pointed out in Ref.[3] the path integral over gauge degrees of freedom 
representing Wilson loops is not of the Feynman type, therefore, it depends explicitly on 
how one "understands" it, i.e. how it is discretized and regularized. We would like to 
emphasize that the regularization procedure applied in Ref. [3] has led to an expression for 
Wilson loops which supports the hypothesis of Maximal Abelian Projection [14]. Accord- 
ing to this hypothesis only Abelian degrees of freedom of non-Abelian gauge fields are 
responsible for confinement. This is to full extent a dynamical hypothesis. It is quite obvi- 
ous that such a dynamical hypothesis cannot be derived only by means of a regularization 
procedure. 
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In order to show that the problem touched in this paper is not of marginal interest 
and to check if path integral expressions that look differently superficially could actually 
compute the same number we evaluate below explicitly the path integrals representing 
Wilson loop for a pure SU{2) gauge field. As has been stated in Ref.[3] for Wilson loops 
C a gauge field along a given curve can be always written as a "pure gauge" and the 
derivation of the path integral representation for Wilson loops can be provided for the 
gauge field taken without loss of generality in the "pure gauge" form. We would like to 
show that for the pure SU (2) gauge field the path integral representation for Wilson loops 
suggested in Ref. [3] fails for a correct description of Wilson loops. Since a pure gauge 
field is equivalent to a zero gauge field Wilson loops should be unity. 

Of course, any correct path integral representation for Wilson loops should lead to 
the same result. The evaluation of Wilson loops within the path integral representation 
Eq.([2.8[) is rather trivial and transparent. Indeed, we have not corrupted the starting 



expression for Wilson loops ( |2.1| ) by any artificial regularization. Thereby, the general 
formula ( |2.8| ) evaluated through the discretization given by Eqs.( |2.7| ) and ( |2.6| ) is com- 
pletely identical to the original expression ( |2.1| ). The former gives a unit value for Wilson 
loops defined for an arbitrary contour C and an irreducible representation J of SU(2): 
Wj(C) = 1. 

Let us focus now on the path integral representation suggested in Ref. [3] 

Wj (C) = f J] DQ(x)e 2U ^cdxMt 3 ^)} y (4 .l) 
J xec 

where all matrices are taken in the irreducible representation J. Following the discretiza- 
tion suggested in Ref. [3] we arrive at the expression 

, „ i , A f s ZiJgfr dx u tr \t 3 A^(x)} , . 

Wj(C) = lim / DQ{x k ) e y jo *m-i»* L ^ v n . (4.2) 

Setting A^x) = d„U(x)W(x)/ig we get 

tf(x) = -d,(Q(x)U(x))(n(x)U(x))l (4.3) 
ig 

By a gauge transformation Q(x)U(x) — > Q(x) we reduce Eq. fl4.lp to the form 

Wj{C) = f J] Dn(x) e 2J f C dxM t3d M x )^(x)} (44) 

J x&C 

For simplicity we consider Wilson loops in the fundamental representation of SU(2), 
Wi/ 2 (C). The result can be generalized to any irreducible representation J. 



For the evaluation of the path integral Eq^^) it is convenient to use a standard s 



parameterization of Wilson loops C [2]: x M — > x^{s), with s G [0, 1] and 2^(0) = x^(l) 
The Wilson loop (|4.4j) reads in the s-parameterization 



1 

w 1/2 (c) = / n on ( s ) ex p / dstr 



0<s<l 



t^nKs) 

as 



(4.5) 
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The discretized form of the path integral Q4.5Q is given by 



As k 



+l,k 



W 1/2 {C) = lim / TT DQ k exp As k+hk tr 

J k=l 

= lim / ft Dn k e tT ^ nk + in ^ = lim /... / DQ n DQ n ^DQ n 

n— >oo J n— >oo J / 
fc=l 

x e tr[t 3 fi n ^_i] e 2Jtr[t 3 O n _ 1 n^_ 2 ] _ e tr[t 3 fi 2 fil] e tr[t 3 ^t 

where fi n+ i = Oi. 

For the subsequent integration we would use the expansion [15] 



-2 



(4.6) 



eM^] =T/aj{z)Xj 



Ant' 



n 



(4.7) 



When the exponent of the l.h.s. is defined for the fundamental representation, the coeffi- 
cients dj{z) are given by [15,16] 



{2 . + 1} 2J 2j+1 (z) 



(4i 



where J 2 j+i(z) are Bessel functions and the index j runs over j = 0, 1/2, 1, 3/2, 2, . . . [16]. 
Recall, that z = 2J and for the fundamental representation J = 1/2 we should set z = 1. 
For the integration over we suggest to use a formula of Ref.[17] modified for our 

case 



(4.9) 



where the coefficients aj(z) are defined by the expansion Eq. ( |4.7| ) . The formula Eq.(O) 
can be derived by using the orthogonality relation for characters [9,10,17] 



2j 



(4.10) 



Integrating over (i — 1, 2, . . . , n) we arrive at the expression 
W 1/2 (C) = lim]T(2j + l) 





n 




Xj 




2j + l 



2j 



}^J2(2j + l)[2J 2+1 (l)} n £e- 



-iirkn 



(4.11) 



fc=0 



where we have denoted 0,(1) = % = e~ m 3 (2j + l)2J 2j '+i(l)- The series over j is 
convergent for any finite n and every term of this series vanishes at n — > oo. This proves 
that the Wilson loop Wi/2(C) vanishes in the limit n — > oo, i.e. Wi/2{C) = 0. 

Thus, the Wilson loop Wi/2{C) for an arbitrary contour C and a pure gauge field 
represented by the path integral derived in Ref. [3] vanishes, instead of being equal to 
unity, Wi/ 2 {C) = 1. This shows that the path integral representation suggested in Ref. [3] 
fails for the correct description of Wilson loops. 
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5 Wilson loop for Z(2) center vortices 

In this Section we evaluate explicitly the path integral ( |4.1|) for Wilson loops pierced 
by a Z(2) center vortex with spatial azimuthal symmetry. Some problems of Z(2) center 
vortices with spatial azimuthal symmetry have been analysed by Diakonov in his recent 
publication [18] for the gauge group SU(2). In this system the main dynamical variable 
is the azimuthal component of the non-Abelian gauge field A^(p) (a = 1, 2, 3) depending 
only on p, the radius in the transversal plane. For a circular Wilson loop in the irreducible 
representation J one gets 

2J + l m tlj 2J+1 sinMp)] 1 ' 



where p(p) = p^jA^(p)A^(p). The gauge coupling constant g is included in the definition 
of the gauge field. For Wilson loops in the fundamental representation J = 1/2 we have 

W 1/2 (p) = cosMp)]. (5.2) 

In the case of Z(2) center vortices with spatial azimuthal symmetry and for the funda- 
mental representation of SU(2) Eq.( |4.2| ) takes the form 

wWp) = lim fr / Dnke tr[t 3 (i2n P /n)n k+1 A^p)nl + t 3 n k+1 ni] = 

k=l J 

= lim / . . . / Dtt n Dn n „ x DQ, n _ 2 . . . DQ X 

n-+oo J J 

x e tr[t 3 (z27rp/n)fi n A (/) (p)^_ 1 + t 3 fi n ^_i] 
x e tr[t 3 (i27rp/n)n n _ 1 A^(p)n] 1 _ 2 + t 3 fi n _ 1 fi] l _ 2 ] 
x e tr[t 3 (z27rp/r2)fi 2 A (/) (p)fi t 1 + t 3 tt 2 {l\] 

x e tr[t 3 ( i 2vrp/n)fi 1 ^(p)^ + t 3 ^] ? (5 3) 

where we have used C Xk+lXk = 2np/n, Q(xk) = &k and fl n+ i = 

For the subsequent evaluation it is convenient to introduce the matrix 

Q{Ap) = (l + i — pA^pj). (5.4) 



n 

In terms of Q^A^) the path integral ( |5.3| ) reads 



^n-l] 



W 1/2 (p) = Jim / ■ ■ ■ / DQnDQ^DQ^ . . . Dtt x e tr[£ 3 O n Q(4^ r 

x e tr[t 3 fi n _ 1 g(A )^_ 2 ] e tr[t 3 fi 2 g(^)fii] e tr[t 3 fi 1 g(A )fit] ; (55) 

The integration over fi^ we carry out with the help of Eq.( |4.9|) taken in the from 

J on k e^l^k+MA^nl + Q(A^Ql_jfn k ] 

2 
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and the orthogonality relation for the group characters. This yields 



Wi/ 2 {p) = lim 



2j + l 



Xj 



(5.7) 



The evaluation of the character Xj[Q n (-^<j>)] f° r n ~ > 00 runs as follows 



( 1 + * — pMp) 



n 



Xj 



sin[(2j + 1)MP)] 
sin[7r/x(p)] 

Substituting Eq.( |5.8| ) in Eq.( |5.7| ) we obtain 



(5.1 



W 1/2 (p)= Hm£[2.Wl)]' 



sin[(2j + 1)MP)] 
sin[7r/i(p)] 



2j 
fc=0 



-mkn 



(5.9) 



The series over j is convergent for any finite n and at n — > 00 every term vanishes. 
This gives W\/2{p) = 0. Thus, we have shown that the the path integral for Wilson 
loops suggested in Ref.[3] gives zero for a field configuration with a Z(2) center vortex, 
W1/2G ) — 0, instead of the correct result W\/2(p) = cos7r/i(p). 

We hope that the examples considered in Sect. 4 and 5 demonstrate that the path 
integral representation for Wilson loops derived in Ref.[3] is erroneous. Nevertheless, in 
Sect. 6 we evaluate explicitly the regularized evolution operator Z Reg (i? 2 , Ri) suggested by 
Diakonov and Petrov for the representation of the Wilson loop in Ref.[3]. We show that 
this regularized evolution operator Z^ eg (R 2 , Ri) has been evaluated incorrectly in Ref. [3]. 
The correct evaluation gives Z Reg (R 2 , R\) = which agrees with our results obtained 
above. 



6 Path integral for the evolution operator Z(i?2,i?i 

As has been suggested in Ref. [3] the functional Z(R 2 ,Ri) defined by (see Eq.(8) of 
Ref. [3]) 



Z(R 2 ,Ri) = J DR{t) exp fir J Tt(iRRr 3 ) 

Ri ^ h 



(6.1) 



where R = dR/dt and T = 1/2, 1, 3/2, ... is the colour isospin quantum number, should 
be regularized by the analogy to an axial-symmetric top. The regularized expression has 
been defined in Eq.(9) of Ref. [3] by 



R2 , t 2 

z Reg (R2, Ri) = J DR(t) exp ii j [- 1 ± (ni + + - J,, nl + m 3 

Ri ^ h 



(6.2) 



where Q a = i Tr(RRr a ) are angular velocities of the top, r a are Pauli matrices a = 1, 2, 3, 
I± and 7|| are the moments of inertia of the top which should be taken to zero. According 
to the prescription of Ref. [3] one should take first the limit 7m — > and then 7^ — >• 0. 
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For the confirmation of the result, given in Eq.(13) of Ref.[3], 
Z Reg (R 2 , Ri) = (2T + 1) D^ T (R 2 R\), 



(6.3) 



where D T (U) is a Wigner rotational matrix in the representation T, the authors of Ref.[3] 
suggested to evaluate the evolution operator ( |6.2| ) explicitly by means of the discretization 
of the path integral over R. The discretized form of the evolution operator Eq. fl6.2p is 
given by Eq.(14) of Ref.[3] and reads^ 

JV 



Z Rcg (R N+ i, Ro) = lim N f ft dR 

N — > OO J n=l 







x exp 



N 



E ( " ^ f( Tr V nn) 2 + (Tr V n r 2 f\ - i i (Tr V n r 3 ) 2 - T (Tr V n r 3 ) 



n=0 



28 



(6.4) 



where R n = R(s n ) with s n = tx+n 8 and Q a = iTr (R n Rl +1 T a ) / 5 is the discretized analogy 
of the angular velocities [3] and V n = i? n i?^ +1 are the relative orientations of the top at 
neighbouring points. The normalization factor M is determined by 



J\f 



N+l 



2ni5 V 2ni5 



(6.5) 



(see Eq.(19) of Ref.[3]). According to the prescription of Ref.[3] one should take the limits 
5^0 and I\\,I± — > but keeping the ratios Ii/5, where (i = ||,-L), much greater than 
unity Ii/5 ^> 1. 

Let us rewrite the exponent of the integrand of Eq. ( |6.4|) in equivalent form 



Z~Reg{RN+l, Ro) — nm 

N — > oo 
5^0 



N+l 



2ni5 V 2iri§ 



N 

U dR r 

n=l 



x exp 



N 



El -*^C&W 



n=0 



i \^ (Tr V n r 3 ) 2 - T (Tr V n r 3 ) 



(6.6) 



Now let us show that if V n is a rotation in the fundamental representation of SU(2), so 



(TrKr a ) 2 = -4 + (TrK) S 



For this aim, first, recall that 



Tr(Kr a ) = -Tr(K f r a ) 



(6.7) 



(6.8) 



Since V n is a rotation matrix in the fundamental representation of SU(2) [19]. By virtue 
of the relation fl6.8[) we can rewrite (Tr V n r a ) 2 as follows 



(Tr Kr a ) 2 = -Tr (Kr a )Tr (V*T a ) = -2 Tr ((V n -~ Tr V n ) " \ Tr K 
= -2 Tr (R n Rl +1 R n+1 Ri) + (Tr K) 2 = -2 Tr 1 + (Tr V n f = -4 + (Tr K) 2 . (6.9) 



x We are using the notations of Ref.[3] 
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By using the relation Eq.( |6.7|) we can recast the r.h.s. of Eq.( |6.6|) into the form 

iii i \ N+l 



lim 

N oo 
5^0 



'4/4)" T ' ex "( l2(Ar + 1) T, 



N 



Yl dR n exp 



n=l 



N 



n=0 



*ll ~ I± 

25 



(TrV^T 3 ) 2 - T(TiV n r 3 ) 



. (6.10) 



Now let us proceed to the evaluation of the integrals over R n (n = 1,2,..., N). For this 
aim it is convenient to rewrite the r.h.s. of Eq. (|6.10 ) in the following form 



^RegC-RiV+l) Ro) — hm 

N -> 00 



JV+l 



x exp 



5^0 

dRfj dRjy—i . . . dR 2 dR\ 
I 



25 



(Tr R N R J N+1 Y + (Tf Rn-^Y + ... + (Tf R 2 R\) 2 + (Tf R 1 R 1 ) 



h-h 



25 



{Tr R n R? n+1 t 3 ) 2 + (Tr R N ^ N r 3 ) 2 + ... + (Tr R 2 R\t 3 ) 2 + (Tr R^^f 



-T 



Tr {R n R} n+1 t 3 ) + Tr (i^i^) + . . . + Tr (R 2 R\t 3 ) + Tr (i2i4r 3 )l . (6.11) 



In the fundamental representation and the parameterization [19] we have 
TrV n = Tr(R n Rl +1 ) = 

fin Pn+1 ( Q-n + ln «n+l + ln+l \ 
= 2 COS COS COS 

2 2 V 2 2 J 

, n . Ai . Ai+l / a « _ 7n "n+l - 7n+l \ 

+2 sm — sm cos = 

2 2 V 2 2 J 

„ / /?n - \ / «n - "n+l \ / 7n ~ 7n+l \ 

= 2 cos cos cos 

V 2 J { 2 J { 2 ) 

[Pn + Pn+l \ ■ fa n - «n+l \ • ( In ~ ln+l \ 

-2 cos ^ j sm ^ j sm ^ j , 

Tr(Kr 3 ) = Tr(i?„4 +1 r 3 ) = 

. Ai Ai+l . / «n + In « n+ i + 7„ + i \ 

= —2 % cos — cos sm 

2 2 V 2 2 J 

. /3 n . (3 n+ i . (a n - 7„ a n+ i - i n+1 \ 

+2 % sm — sm sm = 

2 2 \ 2 2 J 

. [Pn- Pn+l \ ( a n ~ «„+l \ . (in- ln+l 

= — 2i cos cos sm 



[Pn + Pn+l \ . (0L n — a n+1 \ (in- ln+l ' 

-2 z cos sm cos 



(6.12) 
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The Haar measure R n is denned by 



DR n = — ^ sin (3 n d(3 n da n dj n . (6.13) 
87T 

Due to the assumption 1^/8 ^> 1, where (i — ||, J_), the integrals over R n are concentrated 
around unit elements. Expanding Tr (V n ) and Tr (V n T 3 ) around unit elements we get 

TiV n = Tr{R n Rl +1 ) =2-^(/3 n -/5 n+1 ) 2 -i(a n -a n+ i+7„-7„ + i) 2 , 
Tr (V n T 3 ) = Tr (R n Rl +1 T 3 ) = -«(«„- a n+ i + 7„ - 7„+i). (6.14) 
For the subsequent integration it is convenient to make a change of variables 

«n + In — > In, 

C^n 7n //. -,r\ 
~ -> «n- (6-15) 

The Jacobian of this transformation is equal to unity. After this change of variables (|6.14|) 
reads 

TrK = Tr(i? n ^ +1 ) = 2-i(/3 n -/3 n+1 ) 2 -i(7„-7n+i) 2 , 
Tr (V n T 3 ) = Tr{R n R{ +1 T 3 ) = -i( ln - ln+1 ). (6.16) 

Since both TrV n and Tr (V n T 3 ) do not depend on a n , we can integrate out a n . This 
changes only the Haar measure as follows 

DR n = ^f3 n d(3 n d ln . (6.17) 

The integration over (3 n and 7„ we will carry out in the limits — 00 < (3 n < 00 and 

—00 < 7„ < 00. 

Substituting expansions (|6.16|) in the integrand of Eq. ( |6.11| ) we obtain 



ZReg(i?jV+l, i?o) — h m 



/, n^\ N+i a\ N 



oo \2iri5 V 27n<5 y VS 71 ", 
5^0 



x / c/7at / d(3 N p N J d7Ar-i / d/?iv-i /3iv-i • • • J d^ 2 J df3 2 j3 2 J dru J dp 1 (3 X 

-co —00 —00 —00 —00 —00 —00 —00 

x exp (i ^ [(J3 N+1 - (3 N ) 2 + (J3 N - /^-i) 2 + ... + (Jh- A) 2 + (A - A)) 2 ] 
+ i ^ [{iN+i ~ In) 2 + (7jv - 7tv-i) 2 + • • • + (72 - 7i) 2 + (71 - 7o) 2 ] 
- i T [(7jv+i - 7at) + (7w - 1n-i) + • • • + (72 - 7i) + (71 - 7o)] = 



i T (7jv+i - 7b) lim 

JV -»• 00 
5^0 



iV+l / \ N 



h h \ " / 1 



27ri5 V 27ri<5 / \8n 
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oo oo 



oo oo 



oo oo 



x J d'jN J d(3 N (3 N J d-y N ^ J (IPn^Pn-i . . . J d l2 J d[3 2 [3 2 J d 7l J dp!p! 



— CO — oo 

( . I\ , 



—oo — oo 



— oo — oo 



x exp ^ ^ [{p N+1 - (3 N y + ((3 N - /^-i) + ... + ((3 2 - A) 2 + (A - (3oY 
+ i 25 ^ lN+1 ~ + ~ 77V -^ 2 + • • • + (72 - 7i) 2 + (7i - 7o) 2 ] J. 



(6.18) 



The integration over 7„ gives 

00 00 00 00 

J d-f N J c?7 iV _i ... j d>y 2 J c?7i 



—00 —00 



X eXP \26 ^ 7JV+1 ~ + ~ lN '^ 2 + • • • + (72 - 7i) 2 + (7i - 7o) 2 ] J = 



2ixi5 1 



2m5 2 



2?ri5 N - 1 



^ J,, 2^ J,, 3"'^ J,, N ^ J,, N + l 

( h 



2m5 N . h , ^2\ 

eX P I 1 o/Ar (7ah-i - 7o) = 



2(iV+ 1)5 



27Ti5 



N 



N + l 



exp i 



2{N + 1)5 



(7jv+i - 7o) z 



(6.19) 



By taking into account the normalization factor the result of the integration over 7n reads 



(/ ~ v N+l 00 00 00 00 

|m j d lN J d^... J d l2 J dll 



— OO — OO 



OO —OO 



I 



x exp ( % ^ [(7JV+1 - 7iv) 2 + (7iv - 7v-i) 2 + • • • + (72 - 7i) 2 + (71 - 7o) 2 ] J = 



exp 1 



2m(N + 1)5 * V 2(iV + 1)5 



(lN+1 ~ 7o) 



2^ exp i_ ( 7iV+1 - 7 o) 2 , 



(6.20) 



where we have replaced (N + 1) 5 — t 2 — t\ — At. The obtained result is exact. By 
replacing I\\ — > M, 7jv+i — > 70 — > x a and At — > (£& — i a ) we arrive at the expression 
for the Green function, the evolution operator, of a free particle with a mass M given by 
Eq.(2.51) of Ref.[6]. 

Thus, after the integration over 7n the evolution operator Z Rcg (R N+ i, R ) can be 
written in the form 



ZReg(R-N+l, Rq 



2niAt 



exp 1 



2At 



(7AT+1 - 7o) 2 ] e" * T (T^+i - ^ /W ft], (6.21) 



where F[I±, @n+i, Po] is a functional defined by the integrals over j3 n 
F[I±, 0n+i, Po] = 



18 



N+l 



= lim 

5^0 



\ N oo oo 

s) _/«W 



d@N-l Pn- 



dp 2 p 2 / 



x exp \i ^ [(J3 N+1 - (3 N f + (Pn ~ Pn-i? + ... + (&- A) 2 + (A - Pa?] j ■ (6.22) 

Formally we do not need to evaluate the functional F[I±, Pn+x, Pa] explicitly. In fact, the 
functional F[I±, Pn+i, Pa] should be a regular function of variables I±, Pn+i and Pq whose 
absolute value is bound in the limit I± — > 0. Therefore, taking the limit I\\ — > for the 
evolution operator Z^^Rn+i, Rq) defined by Eq.( |6.21| ) we get 



Z(R 2 ,Ri) = lim Z Reg (R N+ i, Rq) = 0. 

in .1 j_ — >0 



(6.23) 



This agrees with our results obtained in Sects. 4 and 5. 

Nevertheless, in spite of this very definite result let us proceed to the explicit evaluation 
of the functional F[I±, Pn+i, Pa] and show that the functional F[I±, P N+1 , Po] vanishes at 
N —>■ oo. It is convenient to rewrite the integrand of Eq. (|6.22j ) in the equivalent form 



F[I±, Pn+i, Pa] 



to f A. fzl 

N ^ oo \2m5) \2n 
5^0 



A' 



8tt 



N 



d d 



d d 



x 



dji dj 2 dj N _x dj 



N 



dPx / dPN 



dp 2 / dp, 



x exp [(P N+1 - p N f + (P N - p N -if + ... + (P 2 - Pxf + (Pt - P ) 2 



3N=3N-l=—=h=jl=0 



+2 j N p N + 2j N „ 1 p N _ 1 + ... + 2 j 2 p 2 + 2 j! p x ] 
After k integrations we get 

oo oo oo oo , 

J dp k J dp k - x ... J dp 2 J dp, exp li^[(p k+1 -p k ) 2 + (p k -p k ^) 2 

— oo — oo — oo — oo ^ 

+ ... + (P 2 - P,) 2 + (Px - po) 2 + 2 3k p k + 2 ]k _x Pk-i + ... + 2j 2 p 2 + 2jx Px]^ = 

eX P \ i n/^u (PO - Pk+l) 2 



(6.24) 



\2m5 1 2m5 2 



<2iti5 k — 1 2ni5 k 



J i 2 V I± 3 



k 



I± k + l 



f.I± a ( k . k 
xexp \i — p k+ x ij-^Jk + 



k-1 



k+l 



■ Jk-i + 



2(k + 1)5 
k-1 k-2 



k+l k k-1 



Jk-2 



+ . . . + 



k 



k-1 k-2 



k-1 k-2 2 1 



k+l 



x exp i 



k 

1 1 

2 ' 2 



72 H 

k-1 3 J k + 1 



k ' /c — 1 " ' 3 ' 2 Jl 



( ■ o\ 2 1 1 o2 1 2 / . 1 . 1x2 



2 3 



1 2 1 

2 ' 3 ' 22 



/3o 2 
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13/ 2. 2 1. 1x2 13 1 a2 14/ 3. 32. 

17l j3 n j2 + r^-^°) + 2-4-3^^~2-5i j4 + 4 J3 + 4-3 J2 

^, 1 5 / 4 43 432 

PO - 7 ' 7 (j5 + 7 J4 + 7 • 7 J3 + 7 " 7 • o 



1 



4 



2 4 

3 2 1 

+ 4"3'2 Jl 

4 3 2 1. 

+ 5 ' 4 ' 3 ' 2 Jl 
k - 1 k-2 



3 2 Ji 3 

2 14 1 



Po) + 

1 \2 

- 5 A) 



4 2 



1 



+ 



fc 



/c - 1 



jfc-2 + • • • + 



2 ' 6 ' 5 2 ^ " 
A; - 1 fc-2 



5 

1 A; 
2 



A- 



A; 

2 . k 
7J2 + - 



5 4 JJ 5 4 3 

fc - 1 . 

3k H ; — Jk-i 



k 



k 



k-2 
k-1 



2 
3 



1 ■ 1 r V 



+- 



A- 



1 



2 fc + 1 A; 2 
By performing N integrations we obtain 



J d(3 N J dp N -!... J d(3 2 J dp 1 expU^[((3 N+1 -p N ) 2 + (p N -p N ^ 1 f 



(6.25) 



+ ... + (02- Pi) 2 + (Pi - A,) 2 + 2j N (3 N + 2 ]N _ 1 (3 N ^ + . . . + 2 j 2 (3 2 + 2 3l fc\ 
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(6.26) 



where we have replaced (AT + 1) 5 = t 2 — t x = At. 

Now we can evaluate the derivatives with respect to ji,j 2 , ■ ■ ■ iJn-IiJn- Due to the 
constraint I±/S ^> 1 we can keep only the leading order contributions in powers of I±/S >> 
1. The result reads 
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Substituting Eq.( |6.27| ) in Eq. ( |6.24| ) we obtain the functional F[I±, Pn+i-, A] : 
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(6.28) 



It is instructive to emphasize that this result can be obtained more easily if we sum up 
like terms in the exponent of the r.h.s. of Eq.( |6.26| ). This yields 
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In order to understand the behaviour of the functional F[I±, 0n+i, 0o] in the limit N oo 
we suggest to evaluate the product 
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at iV 1 by using the £-regularization. In the £-regularization the evaluation of 
II[/3jv + i, /3q] runs the following way 
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Thus, the function H[j3 N+ i, O ] is defined by 
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where e = 2.71828 . . .. Due to the constraint I±/5 ^> 1 the Euler angles 0n+i and 0o are 
less than unity and the ratio in Eq.( |6.32| ) is always positive 
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The functional F[I±, Pn+i, Po] is then defined by 
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Thus F[I±, Pn+i, Po] vanishes in the limit N — > oo. This result retains itself even if we 
change the normalization factor of the evolution operator 
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/ 



vanishes in the limit N — > oo since the Euler angles Pn+i and Pq are smaller compared with 
unity due to the constraint I±/5 ^> 1 [3]. The vanishing of the functional F[I±, Pn+i, Po] 
in the limit N — > oo agrees with our results obtained in Sects. 4 and 5. 
Substituting Eq. (|6.36|) in Eq. ( |6.21|) we obtain 



^Reg(-Rjv+i, Rq) — 0. 



(6.37) 



This leads to the vanishing of the evolution operator Z(R 2 , Ri) given by Eq. (|6.1|) or Eq.(8) 
of Ref.[3], Z(R 2 ,R 1 ) = 0. 

Thus, the evolution operator Z(R 2 ,Ri) suggested in Ref.[3] for the description of 
Wilson loops in terms of path integrals over gauge degrees of freedom is equal to zero 
identically. This agrees with our results obtained in Sects. 4 and 5. As we have shown the 
vanishing of Z(R 2 ,Ri) does not depend on the specific regularization and discretization 
of the path integral. In fact, this is an intrinsic property of the path integral given by 
Eq. (|6.1|) that becomes obvious if the evaluation is carried out correctly. 



7 Evolution operator Z(R2,R\) and shift of energy 
levels of an axial— symmetric top 

In this Section we criticize the analysis of the evolution operator Z(R 2 ,Ri) carried 
out by Diakonov and Petrov via the canonical quantization of the axial-symmetric top 
(see Eq.(12) of Ref. [3]). Below we use the notations of Ref. [3]. 
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The parallel transport operator 

x(t 2 ) 



Pexpfi / ^(x)T a dx M J = P exp fi J A[t) dt 



(7.1; 



a/3 



where = T a dx^/dt is a tangent component of the Yang-Mills field and T a (a = 

1,2,3) are the generators of 577(2) group in the representation T, has been reduced to 
the form 



W aP (t 2 ,t 1 ) = D^(U(t 2 )U\h)), 



(7.2) 



(see Eq.(5) of Ref. [3]) due to the statement [3]: The potential A{t) along a given curve 
can be always written as a "pure gauge" 



d 



A a p{t)=iD T ai {U{t))-D T ^(U\t)) 



(7.3) 



(see Eq.(4) of Ref. [3]). 

By using the parallel transport operator Eq.( [7.2|) the Wilson loop Wt{C) in the rep- 
resentation T has been defined by 



W T (C) = £ W aa (t2, h) = £ D^(C/(t 2 )C/t( tl) ). 



(7.4) 



(see Eq.(25) of Ref. [3]). In terms of the evolution operator Z(R 2 ,Ri) given by Eq.( |6.1D 
(see Eq.(8) of Ref. [3]) the parallel transport operator W a p(t 2) ti) has been recast into the 
form 

W3T(Mi) = // dR x dRa ^2(2T' + l)D^ m (U(t 2 )Rl)DX p (R 1 U\t 1 ))Z(R 2 ,R 1 ), (7.5) 

JJ T',m 

where the index DP means that the parallel transport operator is taken in the Diakonov- 
Petrov (DP) representation. The Wilson loop W^ P (C) in the DP-representation reads 

W^(C)= \\dR,dR 2 £ (2T' + l)D^ m (U(t 2 )Rt)D^ a (R 1 U^t 1 ))Z(R 2 ,R 1 ). (7.6) 

T'm,a 



Of course, if the DP-representation were correct we should get W^ V (C) = Wt{C), where 
W T (C) is determined by Eq.([7~3p. 

The regularized evolution operator Z-& eg (R 2 , Ri) given by Eq.( |6.3| ) (see also Eq.(9) of 
Ref. [3]) can be represented in the form of a sum over possible intermediate states, i.e. 
eigenfunctions of the axial-symmetric top 



Z Reg (R 2 , Pi) = £(2J + 1) Di k (R 2 )D J km (R\) e~ 1 ^ ~ tl} Ejm , 

J,m,k 



(7.7) 



(see Eq.(12) of Ref. [3]), where Ej m are the eigenvalues of the Hamiltonian of the axial- 
symmetric top 



J(J+l)-m 2 {m-Tf 

Ejm = — h 
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(see Eq.(ll) of Ref.[3]). 

As has been stated in Ref.[3]: If we now take to zero Imi — > (first then I±) 
we see that in the sum (12) only the lowest energy intermediate state survives with m = 
J = T. The resulting phase factor from the lowest energy state can be absorbed in the 
normalization factor in eq.(9) since that corresponds to a shift in the energy scale. 

The statement concerning the possibility to absorb the fluctuating factor exp[— i(t 2 — 
ti)T/2ijJ in the normalization of the path integral representing the evolution operator 
is the main one allowing the r.h.s. of Eqs.( |7.5|) and ( fT.6|) to escape from the vanishing in 
the limit I± — > 0. 

In reality such a removal of the fluctuating factor is prohibited since this leads to 
the change of the starting symmetry of the system from SU(2) to U(2). In order to 
make this more transparent we suggest to insert Z-& eg (R 2 ,Ri) of Eq. (|7.7|) into Eq.( [7.6|) 
and to express the Wilson loop W^ P (C) in terms of a sum over possible intermediate 
states, the eigenfunctions of the axial-symmetric top. The main idea of this substitution 
is the following: as the Wilson loop is a physical quantity which can be measured, all 
irrelevant normalization factors should be canceled for the evaluation of it. Therefore, if 
the oscillating factor exp[— i [t 2 — t\) T/2I±] can be really removed by a renormalization 
of something, the Wilson loop should not depend on this factor. 

Substituting Eq . ( |7T7|) in Eq. (|7.6| ) and integrating over Ri and R 2 we obtain the fol- 
lowing expansion for the parallel transport operator in the DP-representation 



W™(t 2 M)=Y: E D^iUit^it^e 



t,)E 



T'm 



(7.9) 



T' m= 



Setting a = (3 and summing over a we get the DP-representation for Wilson loops 



T 



^t DP (C) = £vC P (^i)=E E D^m^it^e-^-^^'m. (7.10) 

a T' m=-T> 

Due to the definition flTJ) the r.h.s. of Eq.f fTloD can be rewritten in the form 

T' 



W™(C) = J2W™(t 2 ,t 1 ) = Yl E WV(C7) e"^ - *0 



(7.11) 



T' m=-T' 



where Wt'(C) is the Wilson loop in the T' representation determined by Eq. ( |7.4|) . The 
relation (|7.11| ) agrees to some extent with our expansion given by Eq.( [7.1|) . 

Following Ref. [3] and taking the limit I\\ — > we obtain m = T. This reduces the 
r.h.s. of Eq.( [7.11[) to the form 



W^{C) = Y J W T ,{C) exp 



-i(t 2 -tx) 



T\V + 1) - T 2 
21, 



Now according to the prescription of Ref. [3] we should take the limit I± 
again Ref. [3] and setting T' = T we arrive at the relation 



(7.12) 
0. Following 



W^(C) = W T (C) exp 



-i(t 2 



T 
2L 



(7.13) 
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Since the average value of the Wilson loop is an observable quantity and any averaging 
over gauge fields does not affect the oscillating factor exp[i(t 2 — ti) T/2I±], one can write 



(W^(C)) = (W T (C)) exp 



- i(h - U 



T 
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According to the Wilson's criterion of confinement Eq.( |1.8|) one should set 



(W T (C)> = e 



-o A 



(7.14) 



(7.15) 



where a and A are a string tension and a minimal area, respectively. 
Substituting Eq. (|7.15|) in Eq . (|7. 14|) one arrives at the expression 
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It seems to be rather obvious that the r.h.s. of Eq.( [7.16|) tends to zero due to a strongly 
oscillating factor, i.e. 
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Really, there is no quantity that can absorb this factor. 

The only possibility to remove the oscillating factor exp[— i(t 2 — ti) T/2I±] is to absorb 
this phase factors in the matrices Ufa) and U\tx) which describe the degrees of freedom 
of the gauge potential A(t) via relation (|7.2|). In this case Eq. ( [7.1 3| ) can given by 



W^(C) = W T (C) exp 



i{t 2 -t x ) 



T 



2I ± 
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2h 



= Y, D L(U(t 2 )U\t 1 ))exp 
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where we have denoted 

U(t 2 ) = U{t 2 )e it ^ T l 21 ^, 

tffa) = U\t 1 )e- itlT l 2lA 



T, D L(U(t 2 )U*(t 1 )), (7.18) 



(7.19) 



However, the matrices U(t 2 ) and W(ti) are now elements of the group U(2) instead of 
SU(2). Thus, the shift of the energy level of the ground state of the axial-symmetric 
top suggested by Diakonov and Petrov in order to remove the oscillating factor changes 
crucially the starting symmetry of the theory from 377(2) to U(2). Since the former is 
not allowed the oscillating factor exp[— i(t 2 — t\) T/2I±\ cannot be removed. As a result 
in the limit I \ — > we obtain 



(W^(C)) 



Reg = lim (W^(C)) = 0. 



(7.20) 



The vanishing of Wilson loops in the DP-representation agrees with our results obtained 
in Sects. 4, 5 and 6 and confirms our claim that this path integral representation of Wilson 
loops is incorrect. 
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8 The non— Abelian Stokes theorem 



The derivation of the area-law falloff promoted great interests in the non-Abelian 
Stokes theorem expressing the exponent of Wilson loops in terms of a surface integral 
over the 2-dimensional surface S with the boundary C = dS [20] 

txVce 19 $c dx n A n{ x ) =trV s e s , (8.1) 

where Vs is the surface ordering operator [20], da^iy) is a 2-dimensional surface element 
in 4-dimensional space-time, x is a current point on the contour C, i.e. x G C, y is a 
point on the surface S, i.e. y G S, and G^ v {y) = d^A u (y) - d u A^(y) - ig[A^{y), A u {y)} 
is the field strength tensor. The procedure for the derivation of the non-Abelian Stokes 
theorem in the form of Eq. (|8.1| ) contains a summation of contributions of closed paths 
around infinitesimal areas and these paths are linked to the reference point x on the 
contour C via parallel transport operators. The existence of closed paths linked to the 
references point x on the contour C is a necessary and a sufficient condition for the 
derivation of the non-Abelian Stokes theorem Eq. ( |8.1| ). 

Due to the absence of closed paths it is rather clear that the path integral represen- 
tation for Wilson loops cannot be applied to the derivation of the non-Abelian Stokes 
theorem. In fact, the evaluation of the path integral over gauge degrees of freedom de- 
mands the decomposition of the closed contour C into a set of infinitesimal segments 
which can be never closed. Let us prove this statement by assuming the converse. Sup- 
pose that by representing the path integral over gauge degrees of freedom in the form of 



the n-dimensional integral (2/7) we have a closed segment. Let the segment C XkXk _ 1 be 
closed and the point x' belong to the segment C XkXk _ 1 , x' G C XkXk _ 1 . By using Eq.(|2.2|) 
we can represent the character x[Ur(C XkXk ^A] by 

x[U?{C XkXk _ x )] = x[U?{C XiX ,)U?{C x , Xi _A] = 

= xMxAUriC^UriC^Mx^)} = 

= d r J DVL r {x')x[Qr{Xi)U r {C XlX >)tll{x')}x[Qr^ = 

= d r [ DQ r (x r )x[U?(C XiX ,)} x[U?{C xlxi _A]. (8.2) 



This transforms a (n — l)-dimensional integral with one closed infinitesimal segment 
into a n-dimensional integral without closed segments. Since finally n tends to infinity 
there is no closed segments for the representation for the path integral in the form of a 
(n — l)-dimensional integral. As this statement is general and valid for any path integral 
representation of Wilson loops, so one can conclude that no further non-Abelian Stokes 
theorem can be derived within any path integral approach to Wilson loops. 



9 Conclusion 

By using well defined properties of group characters we have shown that the path 
integral over gauge degrees of freedom of the Wilson loop which was used in Eq.(2.13) 
of Ref . [1 1] for a lattice evaluation of the average value of Wilson loops can be derived in 
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continuum space-time in non-Abelian gauge theories with the gauge group SU(N). The 
resultant integrand of the path integral contains a phase factor which is not projected 
onto Abelian degrees of freedom of non-Abelian gauge fields and differs substantially 
from the representation given in Ref. [3]. The important point of our representation is 
the summation over all states of the given irreducible representation r of SU(N). For 
example, in 577(2) the phase factor is summed over all values of the colourmagnetic 
quantum number rrij of the irreducible representation j of colour charges. This contradicts 
Eq.(23) of Ref. [3], where only term with the highest value of the colourmagnetic quantum 
number rrij = j are taken into account and the other 2j terms are lost. This loss is caused 
by an artificial regularization procedure applied in Ref. [3] for the definition of the path 
integral over gauge degrees of freedom. 

As has been stated by Diakonov and Petrov in Ref. [3] the path integral over gauge 
degrees of freedom representing Wilson loops is not of the Feynman type, therefore, it 
depends explicitly on how one "understands" it, i.e. how it is discretized and regularized. 
In order to understand the path integral over gauge degrees of freedom Diakonov and 
Petrov [3] suggested a regularization procedure drawing an analogy between gauge degrees 
of freedom and dynamical variables of the axial-symmetric top with moments of inertia 
I± and In. The final expression for the path integral of the Wilson loop has been obtained 
in the limit I±,I\\ — > 0. 

In order to make the incorrectness of this expression more transparent we have eval- 
uated the path integral for specific gauge field configurations (i) a pure gauge field and 
(ii) Z(2) center vortices with spatial azimuthal symmetry. The direct evaluation of path 
integrals representing Wilson loops for these gauge field configurations has given the value 
zero for both cases. These results do not agree with the correct values. 

One can show that Eq. ( |5.9| ) can be generalized for any contour of a Wilson loop in 
577(2) 

Wx/2(C) = fl[lXl{x)e i Bfc dx ^^)] = 

x&C 

(2j + l)W j (C) = 0, (9.1) 

where Wj(C) in the r.h.s. is defined by Eq. fl2.lj ) in terms of the path-ordering operator 
Vc- Further, the result ( |9.1| ) can be extended to any irreducible representation of SU{2). 

This statement we have supported by a direct evaluation of the evolution operator 
^Reg(7?2, Ri) defined by Eq.(14) of Ref. [3], representing the assumption by Diakonov and 
Petrov for Wilson loops in terms of the path integral over gauge degrees of freedom. As we 
have shown in Sect. 6 the regularized evolution operator Z Reg (7? 2 , Ri), evaluated correctly, 
is equal to zero. This agrees with our results obtained in Sects. 4 and 5. In Sect. 7 we have 
shown that the removal of the oscillating factor from the evolution operator suggested in 
Ref. [3] via a shift of energy levels of the axial-symmetric top is prohibited. Such a shift 
of energy levels leads to a change of the starting symmetry of the system from SU(2) to 
£7(2). By virtue of the oscillating factor the Wilson loop vanishes in the limit I± — > 
in agreement with our results in Sects. 4, 5 and 6. 

We hope that the considerations in Sects. 4-7 are more than enough to persuade even 
the most distrustful reader that the path integral representation for the Wilson derived 
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by Diakonov and Petrov by means of special regularization and understanding of the path 
integral over gauge degrees of freedom is erroneous. 

The use of an erroneous path integral representation for Wilson loops in Ref. [7] has 
led to the conclusion that for large distances the average value of Wilson loops shows 
area-law falloff for any irreducible representation r of SU(N). Unfortunately, this result 
is not supported by numerical simulations of lattice QCD [8]. At large distances, colour 
charges with non-zero iV-ality have string tensions of the corresponding fundamental 
representation, whereas colour charges with zero iV-ality are screened by gluons and 
cannot form a string. Therefore, the result obtained in Ref. [7] cannot be considered as 
a new check of confinement in lattice calculations as has been argued by the authors of 
Ref. [7]. 

We would like to accentuate that the problem we have touched in this paper is not of 
marginal interest and a path integral, if derived by means of an unjustified regularization 
procedure, would hardly compute the same physical number as the correct one. We argue 
that no regularization procedure can lead to specific dynamical constraints. In fact, the 
regularization procedure drawing the analogy with the axial-symmetric top has led to 
the result supporting the hypothesis of Maximal Abelian Projection pointed out by 't 
Hooft [14]. Any proof of this to full extent dynamical hypothesis through a regularization 
procedure and through specific understanding of the path integral should have seemed 
dubious and suspicious. 

Finally, we have shown that within any path integral representation for Wilson loops in 
terms of gauge degrees of freedom no non- Abelian Stokes theorem in addition to Eq. flS.ip 
can be derived. Indeed, the Stokes theorem replaces a line integral over a closed contour 
by a surface integral with the closed contour as the boundary of a surface. However, 
approximating the path integral by an n-dimensional integral at n — > oo there are no 
closed paths linking two adjacent points along Wilson loops. Thereby, the line integrals 
over these open paths cannot be replaced by surface integrals. Thus, we argue that any 
non- Abelian Stokes theorem can be derived only within the definition of Wilson loops 
through the path ordering procedure flB.lp . Of course, one can represent the surface- 
ordering operator Vs in Eq. QS.iP in terms of a path integral over gauge degrees of freedom, 
but this should not be a new non-Abelian Stokes theorem in comparison with the old 
one given by Eq. fl8.ip . That is why the claims of Ref. [3-5] concerning new versions of 
the non-Abelian Stokes theorems derived within path integral representations for Wilson 
loops seem unjustified. 
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10 Appendix. Comments on |hep— lat / 0008003] by Di- 
akonov and Petrov 



10.1 None non— Abelian Stokes theorem can be derived within 
path integral representation of the Wilson loop over gauge 
degrees of freedom 

In a set of publications [3,7,15,21,22] Diakonov and Petrov have pointed the possibility 
to derive a non- Abelian Stokes theorem for the Wilson loop represented by a path integral 
over gauge degrees of freedom. This statement contradicts to a theorem that have been 
proved recently by Faber et al. [23] (see also Sect. 8 of this manuscript). The main point 
of this theorem has been proved by using well-defined properties of group characters and 
refers to the well-known result stating that a linear integral over a contour C can be 
replaced by a surface integral (the Stokes theorem) only if the contour C is closed. As has 
been shown in [23] such a necessary condition of the application of the Stokes theorem is 
violated for the Wilson loop represented by a path integral over gauge degrees of freedom. 
Below we repeat our statement. 

In a non-Abelian gauge theory the Wilson loop defined for an irreducible representa- 
tion r of the SU(N) gauge group can be represented in the form of a path integral over 
gauge degrees of freedom as follows [23] (see Eq.(16) of Ref.[23]) 

W r (C) = 4 / II [d r Dn r (x)] X r[U?(C xx )], (10.1) 

where d r is a dimension of the irreducible representation r, Q r (x) is the gauge function 
in the representation r and Xr[U^(C xx )] is the group character defined for the irreducible 
representation r, C xx is a Wilson contour. Then, U^(C XX ) is given by 

U r {C xx ) = V c e l ^c x J^^{z)^ (1Q2) 

where Ay(x) is a gauge field in irreducible representation r of SU(N) gauge group, Vc 
is the operator ordering colour matrices along the path C. Then, the quantity U^{C XX ) 
is defined by [5] 

U?(C XX ) = Q r (x)U r (C xx )Qt(x), (10.3) 

For the matrices U r {C yx ) defined for the open contour C yx linking two points x and y one 
has 

U r (C yx ) = VcJ 9 Sa ~ d ^4 P) ( z ) (10.4) 

and 

U n (C yx ) = n(y)U(C yx )rt(x), (10.5) 

respectively. 
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The evaluation of the path integral over gauge functions Q r (x) can be carried out only 
via the discretization of the path integral ( |10.1| ). Such a discretization can be unambigu- 



ously performed by using well-defined properties of group characters. The discretized 
form reads [23] (see Eq.(15) of Ref.[23]) 

W r (C) = 1 hm J Dn r ( Xn ) XrP^C^J} J M r ( In _,) Xr[U?(C Xn _ lXn _ 2 )} . . . 

x J DQ r (x 2 )xr[U^C X2Xl )] J DQ^XriU^C^J) = 

lim / Dn r (ar ft )Xr[firW^ , r(C r <r„xn-i) n J( a; n-l)] 

x / Dfi r (x n _i) Xr[fy(x n -i)U?(C Xn _ lXn _ a )Stf r (x n - 2 )] . . . 

x [ Dn r (x 2 ) X rin r (x 2 )u r (c X2X M(xi)} f Dn r (x 1 )xrH-(x 1 )u r (c xlWn )nl(x n )], (10.6) 



1 



where C XkXk _ 1 are open infinitesimal segments linking two adjoining points Xk-i and x^ 
and the relations U^}{C XkXk _ l ) = Q r (xk)U r (C XkXk _ 1 )Ql(xk~i) have been used [23]. 

Since infinitesimal contours C XkXk _ 1 are open and the integrations over Q r (xk) are 
independent of the integration over Q r (xk-i) , the necessary condition of the application 
of the Stokes theorem, i.e. a replacement of a linear integral over closed contour by a 
surface integral over an area embraced by the contour, is violated. This implies that none 
non-Abelian Stokes can be derived within the framework of a path integral representation 
of the Wilson in terms of gauge degrees of freedom. 

According to this theorem any non-Abelian Stokes theorem derived by Diakonov and 
Petrov in Refs. [3,7,15,21,22] within path integral representation of the Wilson loop over 
gauge degrees of freedom is very much suspicious and is doomed to be erroneous. 

We argue that in any non-Abelian gauge theory the non-Abelian Stokes theorem 
within the standard definition of the Wilson loop can only be derived via the path- 
ordering operator [20]. This non-Abelian Stokes theorem is unique and no other exists. 

10.2 Path integral representation of the Wilson loop suggested 
by Diakonov and Petrov is erroneous 

In the recent manuscript [15] Diakonov and Petrov have claimed that the path integrals 
over gauge degrees of freedom representing in SU(2) gauge theory Wilson loops defined 
for the pure gauge field and the Z(2) center vortex with spatial azimuthal symmetry have 
been incorrectly calculated in our recent paper [23]. The key point of this claim is the 
incompleteness of the functions Xj[t 3 U], where U is an element of SU(2), U G SU(2), and 
t 3 is the matrix representation of the third generator of SU(2) in the irreducible repre- 
sentation j. Diakonov and Petrov suggested another expansion obeying the completeness 
condition. As has been shown in Sects. 4 and 5 our results for the Wilson loop represented 
by the path integral over gauge degrees of freedom suggested by Diakonov and Petrov are 
retained for the correct expansion. Thus, we conclude that in [15] Diakonov and Petrov 
have not succeeded in refuting our statement that their representation for the Wilson loop 
suggested in Ref. [3] is erroneous and cannot be used for a new check of confinement [7]. 
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10.3 Evolution operator Z Rcg (i? 2 , Ri) is calculated by Diakonov 
and Petrov incorrectly 

First, we would like to cite Diakonov and Petrov (p. 11 of Ref.[15]):"In section 6 of their pa- 
per FITZ0 attempt to compute the regularized evolution operator for the " Wess-Zumino" 
action, following directly our approach. This calculation has been presented in some de- 
tail in original paper [3], however, FITZ seem to be dissatisfied by it and present their 
own. Their final answer (eq.(98) and (112) of ref.fll]), which differs from our, is a result 
of several mistakes. 

First, going from eq.(98) to eq.(87) FITZ use a strange relation, 

ex P ( EH) ^(- 4 )) = ex P {iN(N + 1) y) , (37) 

instead of the correct (and trivial) exp(i2(N + l)I±/5), where I± and 5 are constants and 
iV is the number of pieces in which one divides the contour. 

Second, and more important, both equations in (91) are erroneous, they do not follow 
from eq.(89) from where they are derived. Passing from eq.(89) to eq.(91) one gets: 

Tr(R n Ri +1 ) = 2 - ^[Sa 2 n + 8(5 2 n + 5^ + 25aJ ln cos/3 n ], (38) 

Tr(R n R ] n+1 T 3 ) = i(5a n + 5j n cos/5 n ), 

5a n = a n+i - a n , 5(3 n = (3 n+ i - (3 n , = "f n+1 - j n . 

FITZ have written these formulae without the crucial factor cos/5 n . Because of this 
mistake the subsequent integration over Euler angles a, (3, 7 becomes Gaussian, and the 
evolution operator for the axial top, as computed by FITZ, in fact becomes that a free 
particle, which is definitely wrong. The factors cosf3 n being reinstalled, the derivation 
returns to that of our paper [1]." 

The factor Eq.(37) discussed by Diakonov and Petrov does not effect the final result 
and cancels itself finally. It is rather clear that this "mistake" is nothing more than a 
trivial misprint. Unfortunately, there are some misprints in our paper [23]. For example 
in Eqs.(49) and (74) one has to read Wi/ 2 (C) and Wi/ 2 (p) instead of Wj{C) and Wj(p), 
respectively. In Appendix one should read < a < 2n, < 7 < 2n and < (3 < n 
instead of > a > 2tt, > 7 > 2tt and > (3 > it. 

Before, the explanation of the "Second remark" we would like to attract attention 
of readers to the manner of the expounding of the problem accepted by Diakonov and 
Petrov. 

They write "Second and more important, both equations (91) are erroneous, they do 
not follow from eq.(89) from where they are derived." 

This sentence makes an oritation of the reader that we have made a crucial mistake 
and hardly master the machinery of the expansion of functions in Tailor series. 

After such a successful attack on the psyche of the reader Diakonov and Petrov have 
written completely the same expressions up to the replacement a — > 7 and 7 — > a that 
we have got by expanding the integrand around the saddle point. In addition they have 

2 This is abbreviation from Faber, Ivanov, Troitskaya, Zach used by Diakonov and Petrov 
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added a factor cos (3 n . We are surprised to find such a factor, since in fact there should 
be the factor 



cos 



+ 1 



(10.7) 



This is clearly seen from Eq.(89) of our paper [23]. The variables /3 n +i and (3 n are in- 
dependent. The reason to set (3 n+ i = (3 n would be a mystery of the Diakonov-Petrov 
approach to the evaluation of integrals within a saddle point technique. 

Then, in our paper [23] the factor Eq. (|10.7|) has been replaced by unity. The reason 
of this replacement is in the meaning of the saddle point technique of the evaluation of 
integrals. Indeed, the main point of the saddle point evaluation of integrals is in the 
reduction of integrands to the Gaussian form. In the case of the evolution operator 
ZR eg (i?2j -Ri), as has been claimed by Diakonov and Petrov in Ref.[3], the saddle point 
of the integrand should be at the unit element where a, 7,/? <C 1. That is why the 
expansions should contain the least powers of variables. This imposes the exponent of the 
integrand not to contain the powers of variables higher than 2. That is why the factor 
(|1D.7|) multiplied by (a n -a n+1 ), (a n -a n+1 ) 2 , (7 n -7n+i), (7 n -7 n+ i) 2 should be replaced 
by unity around the saddle point /3 n+ i,/3 n <C 1. 

As we have shown in our paper [23] the evolution operator Z Reg (R 2 , Ri) is equal to 
zero (see Eq.(113) of Ref.[23]) 



10.8) 



This agrees completely with our result W\/ 2 {C) = Wi/ 2 (p) = obtained in this paper. 
Therefore, the claim "The factors cos/3 n being reinstalled, the derivation returns to that 
of our paper [3]." is nothing more than a blef assuming to reanimate the erroneous result. 

Now we would cite again Diakonov and Petrov (see p. 12 of Ref.[15]):" The last ob- 
jection by FITZ is to our alternative (and in fact equivalent) derivation of the evolution 
operator, this time through the standard Feynman representation for the path integral as 
a sum over intermediate states. FITZ quote our result for the evolution operator of an 
axial top with the "Wess-Zumino" term, evolving from its orientation given by a unitary 
matrix R± at time t\ to orientation R 2 at time t 2 : 



ZR Cg (R 2 , Ri) — (2 J + 1 



Djj(R2R\: 



exp 



J 
2h 



(39) 



where I± is a regular moment of inertia, I± — > 0. Apart from a nontrivial dependence 
on the orientation matrices Rl, 2 coming through the Wigner Z)-function, this expansion 
contains a phase factor exp(— i(t 2 — t\) . . .). It is an overall factor independent of the 
external field: it can and should be absorbed into the integration measure to make the 
evolution operator unity for the trivial case R 2 = R\ = 1. Indeed, dividing the time 
interval into N pieces of small length 5, N5 = t 2 — t\, one can write this factor as product, 



exp 



i{t2-h) 



_J_ 
2/7 



A' 

II exp 

fc=i 



iS 



J 

277 



(40) 



where, according to the regularization prescription of ref. [1], S/I± 1 so that each factor 
is close to unity. Each factor can be now absorbed into the integration measure dR(tk) in 
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the functional-integral representation for the evolution operator (39). The fact that the 
factor is complex is irrelevant; moreover, it is typical for the path-integral representation 
of the evolution operators to have a complex measure, see the classical Feynman's book 
[19]. However, FITZ write:". . . a removal of the fluctuating factor is prohibited since this 
leads to the change of the starting symmetry of the system from SU{2) to £7(2)". It 
may seem that FITZ believe that an absorption of a constant factor into the integration 
measure changes the number of variables which one integrates." 

It seems that the decomposition of the finite interval of time into iV pieces t 2 — ti = N 5 
is irrelevant to the problem of the removal of the oscillating factor. Indeed, the expression 
given by the formula (39) does not contain any integration and, therefore, there is no 
measure that can absorb the factor 



exp 



-i(* 2 -*i) 



_J_ 

2L 



(10.9) 



strongly oscillating at Jj_ — > 0. 

In the paper by Faber et al. the crucial contribution of the oscillating factor ( |10.9| ) 
has been demonstrated for the example of the Wilson loop (see Eq.(126) of Ref.[23]) 



Wy F (C) = Wj(C) exp 



-i(t 2 -h) 



_J_ 

2Ti 



10.10) 



where Wf p (C) is the Wilson loop in the Diakonov-Petrov representation and Wj(C) is 
the standard Wilson loop defined via the path-ordering operator. Since the average value 
of the Wilson loop is an observable quantity and any averaging over gauge fields does not 
affect the oscillating factor Eq.( |10.9| ), one can write 

(W? P (C)) = (Wj(C)) exp 7 



2/, 



According to Wilson's criterion of confinement one should set 

(Wj(C)) = e- aA , 

where a and A are a string tension and a minimal area, respectively. 
Substituting Eq. (|10.12| ) in Eq.( |10.11|) one arrives at the expression 

J 



(10.11) 



10.12) 



(W? P (C)) 



e a A exp 



-i(*2-ti) 



21 1 



10.13) 



It seems to be rather obvious that the r.h.s. of Eq. ( |10.13|) tends to zero due to the strongly 
oscillating factor, i.e. 



lim (Wf F (C)) = lim Q e~ aA exp 



~° A lim exp 



i(* 2 -*i) 



_J_ 

2L 



i{t 2 -t t ) 



0. 



J 
2L 



(10.14) 



There is no quantity that can absorb this factor. 

Finally, our explanation concerning the crucial influence of the removal of the oscillat- 
ing factor ( |10.9|) on the symmetry of the gauge fields has been clearly given in our paper 
Ref . [23] . We would not repeat it here and relegate readers to this publication. 
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10.4 Lattice— regularized formula for the Wilson loop suggested 
by Diakonov and Petrov is meaningless 

In this Section we would to make comments on the lattice-regularized formula sug- 
gested by Diakonov and Petrov in Ref.[15]. Formula (50) of the manuscript [15] 

, n 

Wj = AT 1 j dS k exp - Tr{SlU Kk ^S k ^) (50) 

k=l ~ 

is a revised version of the lattice representation of the Wilson loop suggested by Diakonov 
and Petrov in Ref. [7]. However, the parameter z is now not equal to z — 2J, as has been 
stated in Ref. [7], since there is no a priori reason to expect that in the lattice regularization 
z should be the same. [15]. Therefore, now it is suggested to consider this parameter as 
a free parameter z(J) ^ 2.1 . We would like to accentuate that expression (50) for Wilson 
loops according to Diakonov and Petrov is valid for any irreducible representation J of 
gauge group SU(2). 

In order to show that the representation Eq.(50) reproduces the standard Wilson loop 
Wj = 2J X + - Xj{U n ,n-xUn-i,n-2 ■ ■ ■ U liN ) (49) 
Diakonov and Petrov suggested to expand the exponent of (50) according to (34)[|: 



exp zTr[t 3 U] = exp U-iz/2) Tr [e mt ^U] } = E ~ a A z ) Xj 

j 

- 5>(*) f *™bUV)i Mz) = + 1) 

re o 
[20]" [15], 

r N jk 
Wj=N-'J I[dS k j:~a jk (z) E e^D^iSlU^S^), (51) 

k=l jk m k =-j k 



(34) 

j m=~j 

where "the coefficients dj(z) being well-known from the lattice strong-coupling expansion 



Oj [z) = a 



-mj 



2 



(2j + l)-J 2j+1 (z). (52) 



z 



Integrating over the matrices Diakonov and Petrov have arrived at the expression [15] 

Wj = AT 1 J2l b A z )} N xAUN,N^U N ^ N - 2 • • • U 1>N ), (56) 
j 

^ = E(%' + i)[^F. (57) 

j 

where bj(z) = (2/z) J2j+\{z) [15]. Then, Diakonov and Petrov claim that at iV 3> 1 the 
r.h.s. of (56) reduces to the Wilson loop in representation J [15]. 

3 "We make use of the fact that £3 = (— «/2)r3 = (—i/2) expiiirtz) where the last factor is definitely an 
element of SU(2)" [15]. It is obvious that the relation fa = (—i/2) exp^Trt^) is valid only for the matrix 
t3 in the fundamental representation J = 1/2. 
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Now let us adduce some numerical samples demonstrating the erroneous properties of 
the representation of the Wilson loop given by Eq.(56). 

Wilson loop for the fundamental representation J = 1/2, Wi/ 2 (z). For the Wilson 
loop in the fundamental representation J = 1/2 the expansion of the exponential in (50) 
leading to (56) is valid. The result of the numerical calculation of the coefficients in (56) 
and the normalization factor (57) reads 

W 1/2 (z)\ z=3 . 25W3 = 3.027 x 10- 8 + Xi/ 2 (t/ 2 4,23 • • • 

+ 2.264 x Hr 4 Xi(t/ 2 4,23...tA,24) +4.751 x KT 12 X3/2 (^24,23 ••■ U 1>24 ) 

+ ■■■ ~ ^Xl/2(E/24,23---E/l,24)- (10.15) 

The parameter z = 3.25103 is taken from Table 1 of Ref. [15]. Thus, for the fundamental 
representation the expression (56) fits well the standard Wilson loop (49). 
Wilson loop for the adjoint representaion J = 1, Wi(z). First, let us take the 
Diakonov and Petrov point of view and believe that the expression (56) is valid for the 
Wilson loop defined for the adjoint representation. On this way we have to set z = 
4.36765 from Table 1 of Ref. [15] and get the dominant contribution of the character 
Xi (2/24,23 • • • f/1,24) with the prefactor 1/3 according to the normalization of the standard 
Wilson loop (49). The numerical analysis gives 

W^)| 2 =4.36765 = 1.384 X 10 " 9 + 1.795 X 10- 6 Xl/ 2 (t/24,23 ■ ■ ■ U 1>24 ) 

+ 7^ Xl (^24,23 • • • C/-1.24) + 6.366 X 10" 4 X3/2(f/ 2 4,23 • • • 2/l, 24 ) 

+ 1.818 X 10- 10 X2(f/24,23.-.f/l,24)...~^Xl(^24,23---f/l,24). (10.16) 

This means that really the Diakonov-Petrov representation (56) fits well the standard 
Wilson loop (49) defined for the irreducible representation J = 1. 

Wilson loop for the representaion J = 3/2, W^/ 2 {z). Setting z = 5.46564 (see Table 
1 of Ref. [15]) we should get in the r.h.s. of (56) the dominant contribution of the character 

X3/2 (2/24,23 • • • 2/l,24): 

W 3/2 (z)\ z=5A6m4 = 7.145 X lO" 3 + 5.054 X KT 15 Xi/ 2 (f/ 24 ,23 ••• 2/1,24) 

+ 1.430 X 10~ 5 Xl(2/24,23 • • • 2/l,24) + X3/ 2 (2/ 24 ,23 • • • 2/l,24) 

+ 1.089 x 10- 3 x 2 (2/ 24 , 23 • • • 2/1,24) • • • ~ \ X3/2(2/ 24 , 2 3 • • • 2/1,24)- (10.17) 

This fits well the standard Wilson loop (49) defined for the irreducible representation 
J = 3/2. 

Wilson loop for the representaion J = 2, W 2 (z). Setting z = 6.55104 (see Table 1 
of Ref. [15]) we should get in the r.h.s. of (56) the dominant contribution of the character 

X 2 (2/24,23 • • • 2/l,24): 

W 2 (z)\ z=6 . 55104 = 1.073 xlO' 11 + 2.381 xl0- 3 Xi/ 2 (f/ 24 , 2 3---2/i, 2 4) 

+ 2.392 x 10- 22 Xi(2/ 2 4,23 • • • 2/1,24) + 5.096 x lO^M^s • • • 2/1,24) 
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+ ^^X2(f/ 2 4,23 • • • 1/1,24) + 1-5 73 X 10 - 3 X5/ 2 (f/24,23 • • • ^1,2 4 ) • • • 

« ^2(^24,23 •••^1,24). (10.18) 

This fits well the standard Wilson loop (49) defined for the irreducible representation 
J = 2. 

Wilson loop for the representaion J = 5/2, W^/ 2 {z). Setting z = 7.62728 (see Table 
1 of Ref.[15]) we should get in the r.h.s. of (56) the dominant contribution of the character 

X5/ 2 (^/ 2 4,23 • • • f/ 1,24) : 

W 5/2 (z)\ z=7M2728 = 2.051 x 1(T 9 + 3.058 x lO" 7 Xi /2 (£/ 24 , 23 • • • ^1,24) 

+ 3.210 x 10~ 4 X i(U 2 4,23 ■ ■ ■ 1/1,24) + 3.123 x 10 - 38 x 3 / 2 (t/M, 2 3 ■ ■ ■ U 1M ) 

+ 1.201 X 10- 4 X 2 (f/ 2 4,23 • • • 171,24) + X5/ 2 (C/24,23 • • • ^1,24) • • • 

« ^5/2(^24,23 •••^1,24). (10.19) 

This evidences that for J = 5/2 the Diakonov-Petrov representation of the Wilson loop 
(56) fits well the standard Wilson loop (49). 

Wilson loop for the representaion J = 3, W 3 (z). Setting z = 8.69644 (see Table 1 
of Ref. [15]) we should get in the r.h.s. of (56) the dominant contribution of the character 

X3(^/ 2 4,23 • • • £/l,24): 

W 3 (z)\ z=8 m6u = 6.698 x 10- 4 + 1.937 x 10~ 17 Xi/ 2 (£/ 24 , 2 3 ■ ■ ■ U 1<24 ) 

+ 2.681 x 10- 5 Xi(£/ 24 ,23 • • • U 1:24 ) + 2.807 x 10^X3/2(^24,23 • • • U 1M ) 

+ 7.627 x 10- 32 x 2 (f/ 24 , 23 • • • 1/1,24) + 2.216 x 1 0^X5/2 (^24,23 ■ ■ ■ ^1,24) 

+ y^X3(C/24,23---C/l,2 4 ) + ...~ ^X 3 (C/24,23---C/l,24). (10.20) 

Thus, the Diakonov-Petrov representation of the Wilson loop (56) fits well the standard 
Wilson loop (49) defined for the irreducible representation J = 3. 

Now let us explain the real meaning of the Diakonov-Petrov representation for the 
Wilson loop given by (56). This representation is meaningless. 

It is obvious from the following fact. The representation (56) has been obtained via an 
expansion valid only for the fundamental representation J = 1/2. Therefore, (56) should 
be written only as follows 

W 1/2 (z) = M- 1 J2[b j (z)} N X j (U N ,N-iU N - 1 , N - 2 ...U hN ), 

3 

M = E(2j + l)[^Wr- (10-21) 

3 

Without discussing the appearance of the artificial normalization factor Af we argue that 
tuning the parameter z we get the following set of relations 

W l/2 ( Z ) U=3.25103 ~ ^ Xl/2(t/v,Ar-lt/v-l,A^2 • • • U 1>N ) = W 1/2 , 

Wff 2 (z)\ z=i . m7m w i Xi{Un,n-iU N -i,n-2 ■ ■ ■ Ui,n) = W u 
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2=5.46564 ~ 


7 X3/2(Un,N-iUn-1,N-2 ■ 
ZL 


• • C/i,iv) 


= W 3/2 , 


W$(z) 


| z=6.55104 ~ 


1 

7 X2{Un,N-iUn-1,N-2 ■ ■ 



• C/i,jv) = 


w 2 , 


W$(z) 


2=7.62728 ~ 


g X5/2(Un,N-iUn-1,N-2 ■ 


• • C/i,iv) 


= w 5/2 , 


W$(z) 


2=8.69644 ~ 


y Xs(Un,N-iUn-1,N-2 • • 


. u,, N ) = 


w 3 , 



where Wy^(z) is the Wilson loop in the Diakonov-Petrov representation defined for the 
fundamental representation and Wj is the standard Wilson loop for the irreducible rep- 
resentation J = 1/2, 1, 

Since a priori the value of the parameter z is not well defined, so tuning z 
we are able to equate the Wilson loop in the Diakonov— Petrov representation 
defined only for the fundamental representation J = 1/2 to the standard Wilson 
loop defined for any irreducible representation J. This is really a new way to check 
confinement on lattice [7,15]. 



Instead of Acknowledgement 

We conclude that in Ref. [15] Diakonov and Petrov have sold us a conjecture as a 
proof. In a tedious calculation this conjecture turned out to be wrong. A mistake in 
our calculations, unessential for the drawn conclusions, and some misprints were used to 
hide the necessity to redraw the conjecture. New fit parameters were introduced by them 
actually modifying the conjecture. The proof of the new, modified conjecture is again not 
conclusive for J ^ | or reduces to a triviality. We apologize if we fought too hard and 
promise not to recalculate any of Diakonov and Petrov's papers in near future. 
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